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ABSTRACT 


Digital control algorithms have existed for some time in direct 
digital control installations, however, the expense of high resolution 
analog to digital and digital to analog converters makes these 
positional algorithms uneconomic for single loop control. Using 
frequency as a feedback medium, rather than direct A/D and D/A 
conversion, can result in an economical single loop digital 
controller. Employing this concept the traditional position algorithm 
is replaced with a velocity one, in which a voltage controlled 
oscillator and stepper motor replace the A/D and D/A converters 
respectively. Test results from a controller using a PID velocity 
algorithm indicate that frequency feedback can produce a controller 
with performance characteristics similar to those of the conventional 


three mode analog device. 
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Chapter i 


INTRODUCTION 


Between the extremes of the single loop analog controller on the 
one hand and computer control of many loops on the other, lies a field 


of application for an economical single loop digital controller. 


The single loop digital controller could be expected to replace 
the analog device, resulting in better control resolution and drift 
free ait Because the digital controller interfaces easily with 
a digital computer, it would provide an excellent back up for a 
computer operating in the direct digital control environment, or as an 
end device for a supervisory computer system employing setpoint 
control. These obvious advantages of digital control have been 
apparent for some emer? however, the relatively high cost of using 
analog to digital and digital to analog converters with sufficient 
resolution to yield meaningful results has made single loop digital 
controllers uneconomic. An alternate approach in the development of 
the digital controller has been the use of frequency as an intermediate 


medium between the analog and digital states as shown in figure 1. 


WO ae aes i aie ” 
mien ua 


mombgen od treiges ep eee aon ; 
Cie bare fe tidlovc: baie ire sev tly =" 


dite Serine agoeaitst, Satie ins Saoyntn. sith ema F.ryghig nie 
#igett gr daad beet Ciaran ound A 8 eae Kania 

fs he RD yAtmBREVA® Eorthed tabi ti dem teh iid at Sotdiveqo wapigitek S| 
sntagine gutyedget Rabie ‘ae TsO S-\edinh emmgin 8 Se: SRNR SE SE . - 

road Fed Totocs far tpib t akqoemethy etotwte aaalT «eeu | 
gaive Te sabe gt tfovliel We wt rstiocA *opah aby adhd 
sapdetTe-e Wie erottonms guferte' ns scl kon Seca alt 
Det gd accel pitbiede ebiaiy- adit cal Linger Lar kn siete Stes ee ah a 
To. traimpblovah: wit of rigettolge -o7emiolta: oA) tlie mtd. 
sbatecpitns na en ydbeypen® Yo-ain ald pond cbt vara dabigilh 
ie ie or oo a 


7 
_ 
- 
Ml 


- 
: 


7 


digital 
controller 


CONVENTIONAL 


-stepper 
/ motor 
modified 
bts eee digital 
analog controller anajea 
input aaa 
‘Tawi iGanraeireaquency 


FREQUENCY FEEDBACK 


CONTROLLER CONFIGURATION 


FIGURE 1 


Use of the voltage controlled oscillator and the stepper motor as 
analog to frequency and frequency to analog converters respectively 
have eliminated the need for the conventional analog to digital and 


digital to analog devices. 


Problems encountered with transmission of analog signals over 
relatively long distances in a noisy environment can largely be 
eliminated by using frequency feedback. Converting to frequency at 
the end device and using it as the feedback medium provides immunity 
to analog noise and removes the need for extensive filtering to 


produce a" clean " signal. 
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Using the frequency feedback approach Montgomerie et eat produced 
a workable PI temperature controller in 1968. The purpose of this 
research was to extend the scope to include a derivative term, hence 
developing the digital equivalent of the conventional analog three 
mode controller. Several other features were added to enhance the 
controller performance. Instead of providing a control output for 
each of the terms, as was previously done, the components were added 
to produce a composite control signal for each sample period. To 
provide added controller flexibility a sampling rate adjustment has 
been added to allow tuning of the controller sampling rate to the 


process response time. 
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Chapter 2 


THHORY OF OPERATION 


The aim of this project was to implément the digital equivalent 
to the existing analog three mode controller using the frequency 


feedback concept shown in figure 2. 
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FREQUENCY FEEDBACK CONTROLLER 


FIGURE 2 


Consider the stepper motor output shown in figure Shs 
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STEPPER MOTOR CHARACTERISTICS 


FICURE 3 


% 
Providing the step size is sufficiently small the stepper motor may 


be considered an integrator with the following transfer function. 


fs The difference in stepper output and an ideal integrator is shown in 


figure 3. The difference is due to the stepper quantization error. 
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@(t) = Stepper output 
Eat) = frequency input 


to stepper motor 


To produce an output proportional to the measured error the stepper 


motor input frequency must be of the form. 

ad (@) =K,d e(t) e(t) = error (2t2) 

01 1 
dt 
Extending this reasoning, the term corresponding to the integral of 
the error is obtained. 
= ( . 

foo (t) = K, oft) (2.3) 

The derivative term is obtained in a similar fashion. 


2 
£03 (t) = K., d = (2.4) 


Combining the three terms 


fy) (t) = K, 4d oft) + Ky e(t) + K, a” e(t) (2.5) 
at gre 


Using this frequency as the stepper motor input 
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o = K es) 
(c) = KK elt) , KK, feledac , K Kd e(t) (2.6) 
dt 
which is the transfer function of a conventional PID controller. 
Simplifying 
O(c) = Kfe(t) , 1 feledae , T,a et] (2.7) 
Ts dt 
K = KOK 
dae=® K 
i —_—— 
SoKn 
T, = KK 
K 
Since the digital controller is a discrete device, consider the 
discretized equivalent at time t = vee 
=> S= 
aii [ es-1 + 4 2h Vs 9. ae, | (2.8) 
pred At 
At time t = s 
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Assuming no setpoint changes 


ae E 
= = = = A S 
ye Cees [one c.) a at cS cu) en (Ac 4 Ac.) 
Cc 
| | | | 
required proportional integral derivative 
correction term term term 


(2.14) 


The integral and proportional terms can be obtained by a 
straightforward evaluation of the difference between the setpoint and 
feedback signal and between two successive feedback signals 
respectively. Direct evaluation of the derivative term (Ac. 4 - Ac.) 
is difficult; however, it can be evaluated indirectly. It is shown 


in Appendices A and D that the derivative term can be evaluated in 
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both sign and magnitude by measuring the amount of deflection "x" 


in the feedback variable c(t) as indicated below. 


ACg 


c(t) 


ACg_, 


time 


DEFLECTION MEASUREMENT 


FIGURE 4 
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Chapter 3 


IMPLEMENTATION 


To minimize the sampling time, a scheme using four counters as 


shown below was implemented. 
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COUNTER CONFIGURATION 


FIGURE 5 
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Each of the three terms is evaluated individually as follows. 


* 
1 Integral term ae - c.) 


The set point T 5-1 is preset into counter A. The voltage 
controlled oscillator, with c(t) as voltage input, is enabled to 
counter A for a period T. Setting counter A in the count down mode 


produces a residue proportional to (r ae, as shown in figure 6. 


s-1 
The residue in counter A is multiplied by the appropriate gain and 


placed in the accumulator as the integral term. 


: In the interest of minimizing the sampling time, the term 


(x -c_,) will be evaluated instead of ica - c.). The reason 
s-1 s-1 


will become apparent later. 
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INTEGRAL EVALUATION 


FIGURE 6 


2 Proportional term (c4 - c.) 


At time t = t,_, the voltage controlled oscillator with c(t) as 


input is counted up in the initially zeroed counter B for time T. 
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PROPORTIONAL EVALUATION 


FIGURE 7 


After the appropriate delay oe the voltage controlled oscillator is 
counted up in the initially zeroed counter A for period T. Again after 
a delay of &« (at) tine tas t.) the voltage controlled oscillator is 
counted down in both counters A and B to yield the residues shown in 
figure 7. The residue in counter B, which is proportional to 


(cy = ce.) is added to the accumulator as the proportional term. 


+ 
The appropriate choice of «, like T, is dependent on the time 
constant of the controlled system. For added flexibility the terms 


«and T will be variable to allow "tuning" of the controller. 
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3 Derivative term (Ac_, -Ac ) 
. s-1 Ss 


Subtracting one half of the residue in counter B from the 
residue in counter A (see figure 7a), leaves a residue in counter A 
proportional to the magnitude of the deflection of c(t) and hence, as 


shown in Appendix D, proportional to the term (ac. i - Ac). 


residue in A 
after subtraction 
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pulse count 


ts.) ts 
time 


DIFFERENTIAL EVALUATION 


FIGURE 7a 


The counter A residue remaining after subtraction is multiplied by 

the appropriate gain and added to the accumulator as the derivative 
term. Although only one particular case (velocity and acceleration 

negative) 4s shown here, Appendix A demonstrates the consistent 


relationship between the sign of the A residue and (Ac, ~ Ac _) for 
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all possible combinations of velocity and acceleration for c(t). 


To minimize athe the sum of the three control terms stored in 
the accumulator is strobed in parallel into counter C. The evaluation 
of the next sample is started while the contents of counter C is 
clocked to the stepper motor with a proportional gain. Sequence timing 


for a typical sample period is shown in figure 8. 
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Chapter 4 


CONSTRUCTION 


Because most of the construction is a straightforward 
implementation of digital logic, the Tete a the hardware will 
not be discussed but rather only shown on the schematic drawings in 
Appendix C. Only those elements which have a unique construction or 
which directly affect the controller performance will be discussed 


in detail. 


Many design factors interact to affect the controller performance. 
Preliminary considerations suggested that performance Houle be altered 
by such factors as system time constant, sampling rate, quantization 
and gain to name a few. Since this interaction between factors is 
quite complex, mathematical analysis to determine the controller 
settings for the best response is difficult. Ter etereaehe design 
philosophy used in this thesis was to mathematically determine the 
maximum and minimum constraints on the parameters, build the 
controller, then further specify the parameters by evaluating 


preliminary performance tests. 


Several significant components and concepts will now be 


discussed. 


1 Stepper Motor 


Use of the stepper motor as a frequency to voltage transducer 
can be accomplished by coupling the motor output shaft toa 


potentiometer. To enable connection to a standard analog computer for 
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testing purposes, the analog voltage range for both the feedback 

signal into the voltage controlled oscillator and the control voltage 
from the potentiometer attached to the stepper motor was set at 0 to 
+ 10v. Connecting counter C to the stepper motor as shown in figure 9 


results in a digital to analog converter with adjustable gaine 


Counter C is preset from the accumulator at the end of. each 
sample period. The up/down control line for counter C and stepper 
rotation are determined by the sign of the accumulator. Simultaneous 
with the evaluation of the next control signal, ne is enabled to 
both the stepper and counter C and is inhibited when eeenter C 


reaches zero. The amount of stepper movement A is proportional 


t 


to n_ times the original residue in counter C. The stepper used was a 
1.8 yore coupled to a one turn potentiometer which produces a 
quantization error of 1 or + 0.05v. Obviously this error can be 
reduced by employing iit step size or a multiturn potentiometer. 
Performance of this unit will also depend on the stepper input 
frequency (tee) which will affect the rate of output voltage 
change. Testing the stepper motor capability disclosed a maximum 

skip free input frequency of approximately 300 pulses per second. In 
the interest of prolonged stepper life it was desired to reduce this 


frequency, therefore, tee was adjusted during preliminary tests. 


Adequate response speed was obtained for tee = 50 pulses per second. 
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2 Voltage Controlled Oscillator 


Performance of the voltage controlled oscillator is an important 
factor in the performance of the controller. The oscillator frequency 
is counted up in a counter for a specified aperture time T to produce 
a digital result. The amount of resolution, and hence the quantization 
error in this count, is dependent on two factors, namely: the 
oscillator frequency and the aperture time T. Ideally, an infinitely 
high oscillator would produce a good resolution count with a very 
short aperture time. Practically, the linearity of a voltage 
controlled oscillator deteriorates rapidly as the frequency range is 
eee As a compromise, and for reasons of availability, a 


voltage controlled pee ietorun 


with 10 to 10k hertz output was chosen 
(A schematic is shown in Appendix C). The nonlinearity was found to 
be within 0.1%. Having chosen the voltage controlled oscillator, two 


factors remain to be resolved, namely: the aperture and sample times. 


Aperture time 


As shown previously, the discrete nature of the stepper motor 
produces a voltage uncertainty of +0.05v. Because the complete control 
signal consists of three independent components, each component must 
contribute a maximum error of +0.0167v to preserve the output accuracy. 
The errors introduced in each case result from the analog to digital 
conversion and subsequent multiplication. The A/D conversion is 
performed by counting the voltage controlled oscillator frequency, 
therefore,the counter enable time will dictate the quantization error. 


It can be shown (Appendix B) that an aperture time of 0.08 seconds 
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will produce a quantization error of + 0.0125v for each component, 


yielding a worst case composite error of + 0.0375v. 
Sample time 

From the Nyquist criterion the minimum sample rate is defined as 
facie (4.1) 


where fo is the highest frequency component in the sampled signal. An 
accurate evaluation of f. Aspaitficult to obtain, but. cince fo is 
required only to set a minimum bound on fy a rough approximation can 
be made. Because typical response curves are approximately 


Sineergetl, we can consider the sampled curve to be of the form 


c(t) = 2 sin wt (4.2) 
d c(t) = Zw cosw t (4.3) 
at 
d 2) Be) Te (4 4) 


Using the response curve for an analog system with a ten second time 


constant da c(t) ax Was found (figure A-1). 


dt 
dct) es 0.237v/sec = 2w, (4.5) 
dt 
w , = 0-118 rad/sec (4.6) 
f =ew = 0.019 hertz (4.7) 
cayo 
26 


Ba Te 0.038 hertz ae 26 sec (4.8) 


21 


Te ¥ yer a 
: 1 


ot haha be mas. 
ei aa , 7" 
Pa 0 vey PAAAEICY 
nnd come buspar 
~ 


Ae | 
5 7 t.ho la f>« Ths | sa 

(¢.4) wf | aa 

SS 4 

fa, fy | . 

i at? "inane 


Stes) eat, Skt oS pian Be.) 124) 


To allow for the rough approximation used, a sample time of 


approximately two seconds was chosen. 


Making the sample and aperture times variable results in a more 
flexible controller. By accompanying each aperture time increase with 
a corresponding voltage controlled oscillator frequency decrease, the 
final count and hence the resolution, will remain unchanged. This 
feature will allow tuning the controller to a wide range of process 


time constants. The implementation is shown in figure 10. 
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SAMPLE RATE SCHEMATIC 


FIGURE 10 
3 Sequence Generator 


The implementation of the timing diagram shown in figure 8 
required a sequence generator with nine stable states. Because several 
of the interval times resulted from rough approximations, the chances 
of possible alteration, governed by preliminary test results, were 


fairly high. For this reason, an easily adjusted sequence generator 
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was desired. By expanding and modifying the n-flop developed by 
Goldber 16 i 

By a sequencer meeting the above criteria was developed. To 
explain the operation, a partial schematic is shown in figure 11 (A 


complete schematic is shown in Appendix C). 


f —decade counters 
0 
1 
3 o State 1 
4 
5 
6 
7 
8 
9 
oState 2 
7 ine 
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5 | 
6 { | 
Br | | r+ 
G { | i! 
| ' 
( { ', 
a 
1 y to other 
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SEQUENCE GENERATOR SCHEMATIC 


FIGURE iif 


Assuming state 1 is high, feedback through the inverter to all 
NAND gates ensures all other states are zero. The remaining states 
(other than 1) are inverted and used as inputs for the #1 NAND. A zero 


output from the #1 NAND holds state 1 high. State 2 goes high when the 


F bea : 
# g foe ‘ 
7 ps BY a "es 


Ai bP eames. nh 
; - - 15 | a 7 


? 


re in ee 
@islnyos 5 hase 
a aes 


Ph x 


7 ae 
‘ 7 4 } 7 » i { 
Ci TamyOT AS ORE 


counter reaches 06. Feedback from state 2 through the inverter forces 
all other states (other than 2) to zero. This is a stable state and 
will remain until the next desired count is reached. Time duration 
changes are easily made by changing the counter connections to the 
two input NAND gates. Although the generator is extremely flexible, 
it is obvious that in operation one state is always high, therefore, 
in some cases a “dummy" state may be required to produce the desired 
intervals. Another feature is the overlap between intervals. Because 
the next state goes high before the present state goes low, care must 


be exercised when using this sequencer. 
4 Multiplication 


From previous calculations the sample time was chosen as 
approximately two seconds. This slow sample rate allows the required 
multiplications to be carried out serially with very little addition 
to the sample time (see figure 8) and large hardware component 
savings. The programmable divide by n ote is used extensively 
to implement the serial multiplier. A typical configuration is shown 


in figure 12. 
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MULTIPLIER SCHEMATIC 


FIGURE 12 


The value to be multiplied is stored in counter A. The up/down 
count for both the accumulator and counter A is determined by the 
sign of counter A. For example, a positive initial value in counter A 
will set the A counter in a count down mode and the accumulator in a 
count up mode. A zero condition in counter A inhibits both counters 
completing the multiplication. The accumulator count can be shown to 
be 5 times the original count in A. Similar configurations are used 


for each of the three terms in the control equation. Theoretical gain 


calculations for each term are given in Appendix D. 
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5 Sign Convention 


The binary coded decimal number system was used to facilitate 
interfacing to the setpoint input. Positive numbers are represented 
by counting up in the counters yielding the traditional BCD form. 
Negative numbers are obtained by counting down yielding a 1000-x 
(counters are three decade) complement of the original number. A sign 
bit, which is appropriately set by zero transitions, is associated 
with each counter. Each sign is associated with a particular stepper 
motor rotation (positive = clockwise). A negative number in counter C 
is recomplemented on output by counting up to zero with the 


appropriate rotation. The sign convention can be clarified best by an 


example. 
Case 1 Positive 20 
Counter C 
sign bit =1 (positive) 
value = 20 
- sign bit sets the up/down count to down 
and stepper rotation to clockwise. 
-counter C will count down until zero 
detect. Using a gain of 1, stepper will 
rotate 20 steps clockwise. 
Case 2 Negative 20 


Counter © 
sign bit = 0 (negative) 
value = 980 (1000-20) 
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- sign bit sets the up/down count to up and 
stepper motor rotation to counterclockwise. 

- counter C will count up until zero detect. 
Using a gain of Fn stepper motor will 


rotate 20 steps counterclockwise. 
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Chapter 5 


TESTING 


The prime objective of the test series was to evaluate the 
performance of the digital controller in typical control applications. 
Performance was measured by recording the system response to a step 
disturbance in the feedback loop. First and second order lag 
approximations to a process control valve, simulated on the analog 
computer, were used as the controlled systems (see Appendix A). A 
system time constant of ten seconds was chosen as a typical 
representation of the physical system. The test configuration is 


shown physically in figure 13 and schematically in figure 14. 


PHYSICAL TEST CONFIGURATION 


FIGURE 13 
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TEST SCHEMATIC 


FIGURE 14 


Three series of tests, as described below, were conducted to 


evaluate the effect of various parameters. 


Series A (comparison test) 


Using the first order process control valve approximation 
% 
(Appendix E) as the controlled system, the optimized response curves 


for the digital controller were obtained. These curves were compared 


Optimum response, used in the context of this thesis, refers to a 


response that exhibits the best combination of fast rise time and 


minimum overshoot. 
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to the optimized analog controller response for the same system. The 
analog controller used was an analog computer simulation of a three 


mode controller (Appendix E). The results appear in figures 15 - 17. 
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FIGURE 15 
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FIGURE 16 
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Series B (derivative tern) 


The PI control algorithm is optimal for first order sas as 
proven by the series A set of tests (derivative times were approximately 
zero). To evaluate the effect of the derivative component, a second 
order approximation to a process control valve was used (Appendix Ee 
‘Because of the nature of the controller another modification must be 
made. All calculations regarding the derivative component were done 
assuming a fast enough sample rate to ensure that the feat variable 
had a monotonic slope with no inflection points (Appendix A). Using a 
step function as an input it is obvious that this criterion is no 
longer met, thereby invalidating all previous derivative calculations. 
Monotonicity can be restored by expanding the step over the entire 
sample period (figure 1o)e Using this modified step function as an 
input, the effect of derivative gain on the system response was recorded 


in figure 19. 
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Series C (sample time) 


This series of tests demonstrates the well documented mosis 


obtained by varying the Sample time. The results are shown in 


figures 21 - 24, 
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EFFECT OF SAMPLING TOO FAST 


FIGURE 21 
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Chapter 6 


DISCUSSION 


Series A 


Figure 15 is a comparison of the optimum digital and analog 
controller responses to a unit step disturbance. There are two obvious 
differences, namely: the response delay and the steady state 
quantization error of the digital controller, both of which were 
anticipated. Since the curves are virtually identical except for the 
delay, it is interesting to compare the digital and analog gains for 


the components. 
Proportional Gain 


for the digital controller 
K = 0.9 (6.1) 
from Appendix D 
K = 4x” K = 3.6 (6.2) 


for the analog controller 


K = 4.0 (6.3) 


Integral Gain 


for the digital controller 


G, = 0.12 (6.4) 
from Appendix D 

G, Shel) uae 0.056 rep/sec (6.5) 
sok yaaa 
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for the analog controller 


4d = 1 = 0.151 rep/sec (6.6) 


T, 6.6 


Derivative Gain 


for the digital controller 


B = 0.0 (6.7) 
from Appendix D 
% 
T3 = B/p = 0.0 sec (6.8) 
for the analog controller 
T, = 0.1 sec | (6.9) 


Theoretically Ty = 0.0 seconds for optimal control of a first 
order lag process, however, due to the configuration of the analog 
simulation (Appendix E) this was impossible to achieve (required 
infinite gain). Typical analog derivative times range from 0.6 to 
1500 seconds, therefore, Ta = O.1 seconds may be considered zero 


derivative gain. 


The apparently large discrepencies in the proportional and 
integral terms are easily explained by considering the operation of 


the digital controller. The control equation was shown to be 


— “d 
= + = 
Bpoe) the tr eae 3) etoad ihc yea Nee 


T, &t 


(6.10) 


which gives an output correction dependent on gains of the various 
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components. However, using a sample rate several times faster than the 
system response time causes further corrections to be made before the 
system can respond fully to the original. This results in a pseudo 
gain being obtained from the sampling rate. Therefore, the gains 
recorded for the digital controller are somewhat lower than the 
corresponding analog gain as expected. This pseudo gain is a function 
of both the sample time and the system time constant and is difficult 
to compute analytically. Further verification is obtained by 
considering figure 21. The response curve obtained from an excessive 
sampling rate is identical to the response curve which would be 


obtained from an analog controller with significantly higher gains. 


Recall from Appendix D, that to maintain the steady state 
deviation within the desired error bounds, the minimum integral gain 


was limited to G, = 0.25. The integral gain was increased to this 


x 
value and the unit step response recorded (figure 16). We note that 
the transient response shows some deterioration, however, it is now 
ensured that the steady state quantization error is within the limits 
calculated earlier (+ 0.05v). An added advantage of using the higher 


integral gain is the improved setpoint response as described below. 


Another important aspect of the controller operation can be 


demonstrated by considering the control equation. 


T 
ss o + s Tt A as 
o 2 CA K kee oe) At ioe cau d ( Cony e,)| 
T, At 


(6.11) 
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If a setpoint change, rather than a change in the feedback variable 


is produced, this equation reduces to 


6 - Ona K | at (r. - c.)| (6.12) 


This equation results in a response that is less than optimum as shown 
in figure 17. If the response to a setpoint change is critical the 
algorithm must be changed to compensate for this. Mergler et al? have 
developed such an algorithm for a two mode controller. The development 
of an equivalent three mode modified algorithm is beyond the scope of 
this thesis. From the controller timing diagram (figure 8) we note 
that a step disturbance occurring in the to to ty time period will 
also result in the modified control equation (equation 6.12) being 
implemented. Because the occurrence of the disturbance is random, the 
ratio of the to 
of a miss to 12%. Since the controller will be in a physical 
environment with time constants in the order of seconds, a disturbance 
with a rise time of 0.26 seconds (required to cause a miss) will be 


extremely rare. However, if a miss does occur the results will not be 


catastrophic, merely a response that is less than optimun. 


Series B 


The effect of the derivative term on the control of a second 
order system can be readily evaluated from figure 19, which shows 
the response with no derivative term compared with the response for 
T= 0.8 seconds (P=3.0). Because the derivative term obtains its 


a 
effect from the curvature of the sampled signal, the linear input 


to ty period to the total period limits the possibility 
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(modified step) during the first sample period results in no 
improvement during this time. However, from the second sample period 
onward the effect of the derivative term is obvious. Note that the 
addition of the derivative term causes an oscillatory steady state 
which was not present in the first order system. Since the oscillation 
lies within the expected quantization error discussed previously 

(+ 0.05v), this characteristic is not serious. Because of the nature 
of the controller, the derivative term will be effective for 
disturbances with rise times in the order of seconds or longer. For 
disturbances with faster rise times it can be shown, by considering 
the controller operation, that the derivative term will add a 
proportional gain resulting in a response that is less than optimum. 
Since the controller would be expected to operate in a process 
environment where short rise time disturbances are rare or nonexistent, 


the controller performance will not be affected. 


Series C 


This series of graphs demonstrates the effect of sampling rate on 
controller performance. Figure 21 and 22 show the results of sampling 
too fast and too slowly respectively. Figures 23 and 24 show the 


changes in response by decreasing the sample time by a factor of 2. 


The curves are almost self explanatory. A fast sample rate tends 
to over correct for each error because the system does not have time 
to respond to the previous correction before the next sample is taken. 
This produces an oscillatory effect much the same as produced by a high 


gain controller (This pseudo gain is a result of the sample rate as 
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discussed earlier.). Slow sampling allows the system to deviate a 


large amount before a correction is made, yielding much the same type 


of oscillatory result. 


The slower sample rate of figure 24 compared to figure 23, results 


in more delay but still yields an acceptable response. 
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Chapter 7 


CONCLUSION 


As mentioned in the introduction the main objectives of this 
thesis were three-fold: combine the three control Signals into a 
composite control signal; provide added flexibility by employing a 
variable sample rate and probably the most important, verify the curve 
deflection measurement as a derivative control term. From the test 
results it can be concluded that these objectives were achieved with a 


degree of success. 


Throughout the tests, any response deterioration was attributable 
to either the steady state quantization error or the sample rate. 
Because these factors are the common denominator in performance 
limitations, it is beneficial to examine them in more detail. Recall 
that the output quantization error was shown to depend on both the 
step size and the potentiometer used on the output. Therefore, either 
reducing the step size or changing the potentiometer from a single to 
a multi-turn will reduce the error. Changing the potentiometer with 
corresponding changes in gain and stepper input frequency will provide 


the easiest method to improve the steady state error. 


Although the recorded tests were made using a system time constant 
of ten seconds as a typical value, fairly acceptable response was 
obtained for time constants down to four seconds. Operation in 
environments with time constants below this value would have to be 
accomplished by a reduction in the controller sample time. Since it was 


shown that the sample time is a function of the voltage controlled 
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oscillator range, an increase in the output range of the oscillator 
will reduce the sample time. An excellent voltage controlled oscillator 
replacement candidate is the Hybrid Systems Model 610, which would 


replace the original 10k hertz range with a range of up to 100k hertz. 


Making the above modifications coupled with some hardware logic 
simplifications, would produce a cheap high performance digital three 


mode controller for process control applications. 
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Appendix A 


ACCELERATION TERM JUSTIFICATION 


As mentioned in the theoretical portion of the paper, a unique 
relationship exists between the residues of counters A and B and the 
magnitude and direction of the acceleration term. This can be confirmed 
by considering all possible combinations of velocity and acceleration 


as shown in figures A-1 and A-2. 


is Note on figures A-1 and A-2. 
- dashed lines represent negative counter 
residues 
- figure A in each case shows counter residues 
after time period # 6 
- figure B in each case shows counter A 


residue after time period # 7 
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CASE 1 


CASE 2 


Velocity positive (4C,_,-4C,) negative 


Ounter B 
Ounter 


residue 


Velocity positive (4Cs.;-4Cs) positive 


ACCELERATION TERM JUSTIFICATION 


FIGURE A-1 
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CASE 3 


Velo 
city negative (ac, _-4c.) negative 


counter Bg 


counter A 
subtract B/2 
from A 


cece 


pulse count 


negative A 
residue 


CASE 4. Velocity negative (ac._-4c,) positive 


subtract B/2 


pulse count 


ACCELERATION TERM JUSTIFICATION 


FIGURE A-2 
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Summary (Figures A-1 and A-2) 
- add B/2 to A for negative residue in counter A 
- subtract B/2 from A for positive residue in counter A 


Then the sign of the residue remaining in counter A is identical 


to the sign of the term (ac. 4 -ac.). 
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Appendix B 
CALCULATION OF INPUT 


QUANTIZATION ERROR 


The analog to digital conversion is accomplished by counting the 


voltage controlled oscillator frequency for a specified time. 


count = Taw Tt (B.1) 


The counter resolution determines the quantization error as a result 
of converting. Using a voltage controlled oscillator with 10 to 10k 
hertz output for a 0 to 10v input, consider an aperture time of 0.01 


seconds. 


minimum count = 1 (B.2) 
1 = ba 0.01 sec (B.3) 
f = 100 hertz (B.4) 
vco 


Therefore a frequency variation of less than 100 hertz will go 
undetected. This corresponds to an input excursion of 0.iv. Extending 
the above, it is obvious that an aperture time of 0.08 seconds yields 
an input uncertainty of 0.0125v. Because up/down counters are used 
the possibility of either count occurring exists. Therefore, the input 


uncertainty limit for an aperture time of 0.08 seconds is extended to 


if 0,0125v. 
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Appendix C 


SCHEMATIC DRAWINGS 


The following four figures are the complete schematic drawings 


of the components used in the controller construction. 
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CONTROLLER SCHEMATIC 


FIGURE C-1 
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Appendix D 
CALCULATION OF 


THEORETICAL GAINS 


From the theoretical section, the control equation was obtained. 
e@ .-@6 = _ ay 
3 eet K ce c.) + ae (re - c.) P a (Ac 4 a 0c.) | 
i At 


(D.1) 


Proportional Gain K 


The total proportional gain consists of the fixed gains from the 
voltage controlled oscillator and stepper motor and an adjustable gain 


* 
K to allow adjustment of K. 


% 
K=K K_K (D.2) 
K. = stepper gain 


K _ = VCO gain 
vco 


Using a sbefsr /pulse stepper coupled to a one turn potentiometer, one 


pulse yields 1 or 0.05v. 
200 


K_ = 0.05v/pulse (D.3) 


Considering the voltage controlled oscillator with an aperture time of 


0.08 seconds. 


K = 1000 pulse 0.08 sec (D.4) 
Aas volt sec 
= lse/volt (D.5) 
Keg = 80 Pu se/ 
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K = 0.05 volt/pulse ¢ 80 pulse/volt e K (D.6) 
* 
K = 4K (D.7) 


From typical analog data 


aan Peis (D.8) 
Kee E5000 (D.9) 
Yielding 
* 
Kee = 10207 (D.10) 
% 
Kae 255 | (D141) 


% 
To improve the resolution 10K was used as the input switch setting 


(see multiplier construction). 
Integral Gain 


From equation (D.1) the required integral gain is At . From 
th 
/ af ; 
previous calculations (see chapter 4) a sample time of approximately 
two seconds was indicated. For convenience in the sequence generator 


construction, a sample time of 2.15 seconds was used. 


From typical analog controller data 


N 


ay 1500 sec/repeat (D.12) 


i max 


= D.1 
nye 0.6 sec/repeat (D.13) 


These values yield the following gains 
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ae ae at = 2.15 = 1.43010 (D.14) 
1 max 1500 
Ce oe ee SaaS Sys (D315) 
Tt 0.6 
i min 


As shown in the implementation section, the setpoint deviation 
is determined as part of the integral term. Obviously reducing the 
integral term to a near zero value will allow wide deviations from the 
setpoint to occur. Consistent with the expected output accuracy, 
detection and subsequent corrective action for setpoint deviations as 
small as 0.05v is desired. Using an aperture time of 0.08 seconds, 
this deviation yields a counter residue of 4. Since the minimum 
detectable count after multiplication is 1, the minimum gain is limited 


to 0.25. Therefore, the modified gain limits are 


Co ee = eee (D.16) 


= 3.6 (Det?) 


Gs max 


To increase resolution, 10 G, was used as the input switch setting 


(see multiplier construction). 


Derivative Gain 


Pee ¥ 
From equation (D.1) the derivative term is - (Ac. , - Ac.). As 
mentioned in the theoretical section, this term would not be evaluated 


directly but rather the value "x" shown below would be used. 


count 


DEFLECTION MEASUREMENT 


FIGURE D-1 


Appendix A showed the sign relationship between "x" and (Ac. _, = Ac.). 


The magnitude relationship remains to be found. 


s-l cas Ac. (D.18) 


where B is to be determined 


For convenience in evaluation set ue = 1 second 


then 


p*|x | = ny eee (D.19) 
At 
% 
p. ty°F (D.20) 


Consider a constant acceleration curve of the form 
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e(t) = 4 at” 
a 


(D.21) 


count 


DERIVATIVE TERM EVALUATION 


FIGURE D-2 


Determine the difference in slopes. 


=, d c(t) (De22) 
At cea pene cep |e 
From (D.21) 
d c(t) = at (D.23) 
dt 
Ac = Con te a 
s-1 Sd [oe el art ts) (D.24) 
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Determine |x| 


S-1 


|x| = P - c(t) 
ra tot 
= s s-i 
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Substituting equations (D.24) and (D.29) into (D.19) 
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(D.31) 


| i Sen 
‘ : - ‘a “4 . 
ii le wel *. 


$ | . a 
c Bao Be, 
(ds) (raul apna (Ge 
‘ U - | ; i 
' : oly 
; ( : f a 
te 7 f _ ‘ aa Me 
foeon) ™ ~ fof 
eS ; f e r att i| ; ; a \ i‘ 
ny ; fis : 'y) nn 7 
i dee PP a ela 


a 1’ i ¥ : | 
(eta) ogat (out) poe (Git) edottsuge gute hia 
é i ‘ 4 4 

| : 


i 


2 
rk ¢ 
' 


> i, 


Using a sample time of 2.15 seconds 


* 
f = 3.72 
From typical analog controller data 


Ae ee 1500 seconds 


ik = 0.6 seconds 


dad min 


From equation (D.20) 


B= Togs 
B nay = 1500 3.72 = 5580 
[9 MUSE Cys ees 
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(D.35) 


(D.36) 
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Appendix E 


ANALOG COMPUTER DIAGRAMS 


The following two figures represent the analog computer 
simulations for both the analog controller and the system (control 


valve) used in this thesis. 
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SERIES A- FIRST ORDER LAG 


H(s) _ G 
Q(s) 1+Ts 


SERIES B - SECOND ORDER LAG 
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CONTROL VALVE SIMULATIONS 


FIGURE E-2 


67 


| 

_ ba : 

Ce Rad 
" DAS azar 


‘ - 


; wa fy tw ae 7 Vw 
7 amen 
Turd ye dts es 7 


mf ee 
a? i 
» 


a ihe. 


a ines 


ae aa _ ce ah 


an th 
Pt via 


Ob vay ikea 
ia 


Te» <-% 


wun on a ae.  , . : P “a 7 
Th ie, Oe, ne ee arr | 


